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Diffuse  reflected  light  is  used  to  characterize  optical 
properties  of  two-layer  turbid  media.  The  source  is  mod¬ 
ulated  in  gigahertz  frequencies.  The  phase  and  amplitude 
of  analytical  reflectance  is  fitted  to  the  data  generated  by 
adding  2.5%  noise  in  amplitude  and  1°  noise  in  phase  to 
the  solution  of  diffusion  approximation  in  gigahertz  fre¬ 
quency  regime.  The  extracted  optical  properties  are  in 
good  agreement  with  the  real  values. 
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I.  Introduction 

There  are  many  techniques  to  image  living  tissues  non- 
invasively  such  as  anatomical  (MRI,  X-ray...)  and  bio¬ 
chemical  (PET,  fMRI,  optical  methods).  Both  imaging 
modalities  are  used  in  a  complimentary  way  to  distinguish 
an  abnormal  tissue  from  the  rest  of  the  healty  tissues. 
While,  x-ray  imaging  employs  ionizing  radiation,  MRI  is 
very  expensive  and  not  portable.  PET  requires  the  in¬ 
jection  of  a  radioactive  substance  into  the  blood  stream 
and  emission  is  mapped;  hence  it  is  invasive  and  expensive 
system  as  well.  Therefore,  an  inexpensive,  portable  appa¬ 
ratus,  which  can  provide  both  structural  and  dynamical 
changes,  is  a  very  attractive  feature.  Near-infrared  (NIR) 
light  techniques  fulfill  these  requirements. 

In  NIR  region,  tissue  optical  properties  (absorption  and 
scattering)  are  being  mapped.  Light  absorbing  and  scat¬ 
tering  structures  carry  information  about  tissue  function 
and  structure.  Therefore,  quantification  of  optical  prop¬ 
erty  changes  gives  information  about  disease-related  func¬ 
tional  and  structural  changes. 

The  probing  depth  in  the  tissue  depends  on  the  modula¬ 
tion  frequency  of  the  source.  Because  of  the  attenuation, 
at  higher  modulation  frequencies,  the  penetration  will  be 
at  the  surface,  or  skin  layers.  It  is  reported  that  for  modu¬ 
lation  frequencies  between  several  hundred  megahertz  and 
a  few  gigahertz,  the  properties  of  muscle  tissue  of  up  to 
4-8  mm  deep  can  be  analyzed  [1].  It  is  suggested  that 
the  technique  might  be  of  use  for  evaluating  the  depth  of 
necrosis  and  blood  volume  fraction  in  deep  burns.  The 
non-invasive  determination  of  the  depth  of  severe  burns 
has  been  used  for  clinical  applications. 

Tromberg  et  al.  used  an  apparatus  which  goes  up  to  1.5 
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Fig.  1.  Geometry  of  two-layer  model 

GHz  modulation  frequencies  to  determine  optical  prop¬ 
erties  of  two-layer  problem  and  they  suggest  that  going 
even  higher  than  2  GHz  may  improve  the  accuracy  of  ex¬ 
tracting  all  five  free  parameters  in  a  two-layer  system  (ab¬ 
sorption  and  scattering  coefficients  for  each  layer  +  thick¬ 
ness  of  the  layer)  [2].  This  is  due  to  the  increased  scat¬ 
tering  sensitivity  and  better  spatial  localization  of  high- 
frequency  diffuse  photon  density  waves  (DPDW’s).  The 
advantages  of  gigahertz  DPDW’s  are  that  they  are  more 
sensitive  to  boundaries,  and  therefore  in  reconstruction, 
the  artifacts  coming  from  the  mismatch  at  the  boundaries 
will  be  eliminated  [3]. 

In  section  II  we  show  the  analytical  expression  for  two- 
layer  reflectance.  In  section  III  we  obtain  the  optical  pa¬ 
rameters  by  fitting  the  reflectance  of  theory  to  the  exper¬ 
iment. 


II.  Theory  of  Two-Layer  Reflectance 


The  investigated  turbid  medium  is  shown  in  Fig-1. 

The  thickness  of  the  first  layer  is  L,  and  optical  pa¬ 
rameters  are  (/iao,/i'so),  for  absorption  and  scattering  co¬ 
efficients,  respectively.  Second  layer  is  semi-infinite  with 
optical  parameters  (tiai,  /4i)-  It  is  assumed  that  the  two 
layers  have  the  same  refraction  index  n  =  1.33.  The  dif¬ 
fusion  equation  in  the  gigahertz  regime  is  [4]: 


1 

DV2U(r,t)  -  n*U(r,t)  =  (1  +  3 D»a)--U(r,t) 

c  ot 

3 D  d2U(r,t) 
q+~^  dt 2 


(1) 


where  U(r,t)  =  U  (r)exp(iut)  is  the  average  intensity  and 
q  is  the  source  term,  c  is  the  speed  of  light  and  D  is  the  dif¬ 
fusion  coefficient.  When  a  point  light  source  is  modulated 
with  a  frequency  co,  U(r,t)  represents  gigahertz  DPDW’s 
with  a  wave  vector,  k0  =  (—jia/D  +  f|g(  1  +  3 jj,aD)  + 
3  #n2)1/2. 

At  any  constant  z-plane,  U(r)  can  be  expressed  in 
angular  spectrum  representation  as  a  superposition  of 
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plane  waves  with  amplitude  A(K)  and  wave  vector  k  = 
(■ K,q),\k\  =  k0  [5]: 

U(R,  z)=  J  A(K)exp(iK  •  R  +  iq\z\)dK  (2) 

where  R  =  (x,y),  K  =  ( Kx,Ky ),  and  K2  +  q2  =  k2. 

Because  it  is  simple  and  physically  illustrative,  we  solve 
the  two-layer  problem  by  successively  adding  multiple  re¬ 
flections  and  transmissions  at  the  interfaces  [5,6].  Consid¬ 
ering  all  multiple  reflections  from  the  interfaces,  the  total 
reflection  coefficient  at  z  =  0  is  : 


UR  =  UITT , 


10 


+  U. 


R10exp(2iq1L)R12Ti 


10 


zl 


J  zl 


1  -  Ri0exp(2iqiL)Ri2 
R12exp(2iq1L)T10 
1  -  Rwexp{2iq1L)R12 


(3) 


where  Z/®"^  are  the  incident  fields  at  zl  and  z2  respec¬ 
tively,  Rij  and  are  the  reflection  and  transmission  co¬ 
efficients  on  going  from  medium  i  onto  medium  j  and  q  is 
as  before: 

q{K)  =  \jhl-\K\*  (4) 


III.  Fitting  Results  and  Conclusion 
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Experimental  data  is  obtained  by  simulating  the  for¬ 
ward  solution  of  diffusion  equation  with  an  addition  of 
amplitude  noise  of  2.5%  and  phase  noise  of  1°,  because 
instrumental  noises  are  in  this  range  [2].  The  optical  pa¬ 
rameters  are  taken  as  the  same  as  Tromberg  et  al.  for 
comparison  reasons  [2].  The  simulated  data  are  fitted  to 
the  analytical  result  for  reflection  coefficient,  Eq-4.  To 
perform  fitting,  Levemberg-Macquart  Matlab  fitting  al¬ 
gorithm  is  used  to  extract  the  optical  parameters  of  sec¬ 
ond  layer,  because  in  most  experimental  setups  first-layer 
properties  are  already  known.  Table-1,2  shows  the  fitting 
results  for  modulation  frequencies  of  ui  =  1, 2 GHz.  The 
thickness  of  the  top  layer  is  taken  as  L  =  5 mm.  As  seen 
from  the  table  that  the  fitted  values  are  in  good  agreement 
with  the  real  values. 


TABLE  I 


u  =  1  GHz 


Coefficients 

Med-1 

Med-2 

Fitted  Med-2 

fJ*a 

0.05 

0.056 

0.0515 

P-'s 

3.1 

6.2 

6.226 

TABLE  II 


u = 2 GHz 


Coefficients 

Med-1 

Med-2 

Fitted  Med-2 

fJ*a 

0.05 

0.056 

0.0637 

P-'s 

3.1 

6.2 

6.149 

